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ABSTRACT 

We present results of the instability analysis of the post-main sequence massive star models 
against radial and nonradial pulsations. We confirm that both p- and g-modes can be excited 
by the /(--mechanism acting in the metal opacity bump. However, as opposed to the previous 
claims, we find that an intermediate convective zone (ICZ) related to the hydrogen burning 
shell is not necessary for excitation of g-modes. These modes can be reflected at a minimum 
of the Brunt-Vaisala frequency, located at the top of the chemical composition gradient region 
surrounding the radiative helium core. This minimum is associated with the change of actual 
temperature gradient from the adiabatic value in the semiconvective zone to the radiative 
value above it. Thus, the existence of pulsations at this evolutionary stage does not prove 
the existence of the convective zone but only some reflective layer. Finally, we show that no 
regular patterns can be expected in oscillation spectra of blue supergiant pulsators but there is 
a prospect for identification of the mode degree, (, from multicolour photometry. 
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1 INTRODUCTION 

The properties of pulsations of the main sequence B-type pulsators 
such as fi Cephei stars and Slowly Pulsating B-type stars (SPB) 
have been studied extensively for years, but very little is known 
about pulsational properties of stars after Terminal Age Main Se- 
quence (TAMS). Moreover, until recently no pulsating post-SPB 
star has been found. Using the data from the MOST satellite for 
the blue supergiant HD 163899 (B2 Ib/II, Klare & Neckel 1977, 
Schmidt & Carruthers 1996), Saio et al. (2006) found 48 frequen- 
cies with the values less than 3 c/d and maximum amplitudes of a 
few millimagnitudes. The authors attributed these light variations 
to the g- and p-mode pulsations. HD 163899 is a less luminous star 
than a Cyg variables (van Leeuwen et al. 1998) and the first of 
previously unknown type of variable, i.e., Slowly Pulsating B-type 
supergiants, (SPBsg), as termed by Saio et al. (2006). 

According to previous calculations (Pamyatnykh 1999; Pamy- 
atnykh & Ziomek 2007) the instability domain of SPB stars was 
perfectly confined to the main sequence and the instability band 
of/? Cephei stars was slightly extended behind TAMS. Saio et al. 
(2006) examined pulsational instability of models in the mass range 
of 7 < M/M Q ^ 20. The most important result of their analysis was 
to find the large g-mode instability domain in post-main sequence 
models. In the hotter part of this instability region, p-modes are also 
excited. It is in agreement with stellar parameters estimated for HD 
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163899 which locate the star in this instability band. Moreover, ac- 
cording to computations of Saio et al. (2006), the frequency range 
of HD 163899 covers the g- and p-mode pulsations. 

Before the discovery of pulsational modes in HD 163899 it 
had been believed that there is no possibility for g-modes to be ex- 
cited in B-type supergiants. It is due to a strong radiative damping 
which occurs in a dense radiative helium core due to the very large 
value of the Brunt- Vaisala frequency. Saio et al. (2006) have ex- 
plained excitation of g-modes in such stars owing to a partial reflec- 
tion of some pulsational modes at an intermediate convective zone 
(ICZ) related to the hydrogen burning shell. Thus, ICZ can prevent 
some modes from entering the radiative core and /e-mechanism in 
the superficial layers might be effective. This result has been con- 
firmed by Godart et al. (2009) who studied also the effects of mass 
loss and overshooting from the convective core during the main se- 
quence phases. Lebreton et al. (2009) examined the effect of the 
adopted convective instability criterion. All these authors claimed 
that formation and properties of ICZ was crucial for a presence of 
unstable modes in SPBsg stars. However, in this paper, we show 
that pulsational modes in the B-type supergiant models can be re- 
flected at the chemical composition gradient formed above the ra- 
diative helium core and a presence of ICZ in the hydrogen burning 
shell is not conditio sine qua non. 

Here, we are not going to interpret the oscillation spectrum of 
HD 1643899 because the effective temperature, luminosity and ro- 
tational velocity of the star are very poorly determined or unknown 
and we postpone this task for the future when better observations 
will be collected. 

The structure of the paper is as follows. In Section 2, we 
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present the propagation diagram for a representative model and do- 
mains of pulsational instability for SPBsg stars on the HR diagram. 
Properties of instability parameter and kinetic energy of modes are 
described in Section 3. Section 4 is devoted to the photometric di- 
agnostic diagrams as tools for identification of the SPBsg modes. 
The last section contains Conclusions. 



2 EXCITATION OF OSCILLATION MODES IN B-TYPE 
SUPERGIANTS 

Regions of a star where pulsation modes can propagate in radial 
directions are strictly defined by the values of the Lamb (L?) and 
Brunt- Vaisala (A' 2 ) frequencies, which are give by 
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respectively. In the above formulae, I is the mode degree, c s is the 
sound speed, V is the temperature gradient and V flrf its adiabatic 
value, V^ is the mean molecular weight gradient, and (3 is a ra- 
tio of the gas pressure to the total pressure. Other symbols have 
their usual meanings. Wave propagates radially with the frequency 
a) only if a) 2 > L 2 ,N 2 (p-modes) or a> 2 < L 2 ,N 2 (g-modes) and 
they are evanescent in the other cases. This results from the disper- 
sion relation which relates the pulsational freqeuncy with the radial 
wavenumber. As can be easily seen from the above formulae, the 
Lamb frequency depends mostly on the mode degree and the dis- 
tance from the center, whereas the Brunt- Vaisala frequencies on the 
stellar structure details; in particular, on a presence of the chemical 
composition gradient and convection. 

Here, we will consider B-type supergiant models in the hydro- 
gen burning phases. Such star is made up of the radiative helium 
core, the hydrogen burning shell and extended envelope. The other 
possibility is that these objects have already entered the core helium 
burning phase but it is not essential for the point. 

The evolutionary models were computed using the Warsaw- 
New Jersey code. We adopted both the OP (Seaton 2005) and 
OPAL (Rogers & Iglesias 1996) opacity tables. The AGSS09 chem- 
ical mixture (Asplund et al. 2009) was used with the hydrogen and 
metal abundance by mass of X = 0.7 and Z = 0.02, respectively. 
In the Warsaw-New Jersey code the Ledoux criterion for a con- 
vective instability is applied, which states that convection appears 
if N 2 < 0. In unstable layers the adiabatic temperature gradient is 
assumed but the element abundance is frozen. Justification for this 
procedure was given by Dziembowski (1977). All our computa- 
tions were done in the framework of the zero-rotation approxima- 
tion, i.e., all effects of rotation on stellar structure and pulsations 
were neglected. Moreover, no overshooting from the convective 
core during the main sequence phase was allowed. 

An example of the propagation diagram for a model with a 
mass of 16M G is shown in Fig. 1. The model was computed with 
the OP opacities and has effective temperature of log r cff = 4.340 
and luminosity of log L/L Q = 4.730. Due to high mass concentra- 
tion in the radiative core, the value of N is very large in the center. 
The Lamb frequency is plotted for the mode degree I =1, 2. Both 
frequencies are in units of [c/d] . The radiative helium core is very 
dense and of a small size; its border is at the temperature of about 
log T = 7.6 which corresponds to about 4% of the stellar radius. 
Just above the core, we have the chemical gradient zone, whereas 
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Figure 1. Runs of the Brunt-Vaisala frequency, N (solid line), and the Lamb 
frequencies, Lg, for I = 1,2 (dashed and dotted lines respectively) in [c/d] 
for a star in the supergiant phase with a mass of M = 16M ,logr c ff = 4.340 
and log L/Lq = 4.730. The model was computed with the OP opacities 
and AGSS09 solar mixture was adopted. The frequency ranges of unstable 
modes with t = 0, 1,2 are marked as dashed stripes. 
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Figure 2. The hydrogen profile and runs of the temperature gradients in the 
model considered in Fig. 1. 



convection connected with the hydrogen shell burning has not been 
developed because N 2 > in this region. This is due to the fact 
that in regions of nuclear burning the values of ji increases with 
increasing depth and the gradient V^ contributes positively to N 2 . 
Dashed strips show the frequency ranges of unstable modes with 
I = 0, 1, 2. It is worth to note that in this model the radial funda- 
mental mode is unstable whereas the higher overtones are stable. 
A model with similar parameters but computed with the OPAL 
data has a similar overall run of /V. The difference is that the OPAL 
model reaches the lower minimum of N around log T = 7.1 and 
has thinner convective layer associated with the Z-bump. To give 
a better insight into the internal structure of the considered model, 
in Fig. 2, we plotted the hydrogen profile, X, and three temperature 
gradients: actual (V), adiabatic (Vaj) and radiative (V ra( i). As we can 
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Figure 3. Instability parameter, r\, as a function of frequency for I = 0, 1, 2 for the post-MS (left panel) and MS model (right panel) with the mass of 16M© 
and log r c ff = 4.340 and 4.445, respectively. The models were computed with the OP tables and AGSS09 mixture. 
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Figure 4. The HR diagram with instability domains, marked as thick lines, for the radial (upper panels) and nonradial (lower panels) modes of models with 
masses of 7 - 20M Q . The left panels show models with unstable modes computed with the OP opacities whereas the right panels - with the OPAL data. The 
AGSS09 chemical mixture, hydrogen and metal abundance of X = 0.7 and Z = 0.02, respectively, were assumed. 
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see, almost in the entire interior, we have V = V rad , i.e., the whole 
energy is transported by radiation. Only in the Z opacity bump at 
log T ~ 5.2, convection is quite effective and carries up to 25 % of 
the flux energy. Our evolutionary models have been computed with 
the value of the mixing-length parameter, q-mlt, equal to 1.8. The 
choice of this value influences the quantitative results on the pulsa- 
tional instability, but qualitatively the results remain very similar, as 
we tested using the models computed with o-mlt = 0.5. In the chem- 
ical gradient zone (semiconvective zone), we have V ra(i > V ad but it 
does not mean that some energy must be carried by convection. As 
it was discussed by Dziembowski (1977), the required part of en- 
ergy can be carried by short-wave nonadiabatic oscillations. Mixing 
of chemical elements is not associated with this process. Therefore 
in the stellar model calculations we use the simplest description of 
the semiconvective zone structure, i.e., we assume V = V ad and do 
not mix the matter in this zone. Convection here is stabilized by the 
hydrogen abundance gradient. Therefore, the Brunt- Vaisala (N 2 ) 
frequency is greater than zero there. Note also that a local mini- 
mum of the radiative gradient at logT of about 7.6 is connected 
with inner boundary of the hydrogen burning shell, where luminos- 
ity quickly increases outwards. Such a feature does not influence 
dynamical properties of the interior which we discuss in the paper. 
In the next step, we have checked instability of supergiant 
models against radial and nonradial pulsation. To this end, we used 
the non-adiabatic code of Dziembowski (1971, 1977). As one can 
see in Fig. 1, two widely separated regions of unstable modes are 
present. The one occurs at low frequencies (0.2-0.3 c/d) and corre- 
sponds to gravity modes with very high radial orders, n a 250-350. 
The second one at frequencies 1.8-2 c/d is related with lower order 
g-modes (n a 30 - 50) which behave like p-modes only in very 
outer layers and with a radial fundamental mode. These results are 
presented also in the left panel of Fig. 3 where we show the insta- 
bility parameter, 77, as a function of the mode frequency for modes 
with I = 0, 1 and 2. For a comparison, in the right panel of Fig. 3, 
we put a typical run of r]{v) for a main sequence model with the 
same mass and logr cff = 4.445. Let us recall that 77 measures the 
net energy gained by a mode during one pulsational cycle and is 
defined as 

W 
1=~, ■ O) 



I dr 



dr 



where W is the global work integral. If 77 > 0, driving overcomes 
damping and a mode is unstable. As we can see the general depen- 
dencies of rj(v) for the post-main sequence (post-MS) model and 
main sequence (MS) model are similar, i.e., there are two global 
maxima of 77. The first one is related to the high-order g-modes and 
occurs at v a 0.9 and 0.25 c/d for the MS and post-MS model, re- 
spectively. The second global maximum corresponds to low order 
p/g-modes and occurs atv = 6.2 and 1.8-2.0 c/d for the MS and 
post-MS model, respectively. The difference is that for the super- 
giant model, within the low frequency global maximum there is an 
additional pattern of many local maxima and minima, as we have 
already mentioned, in the post-MS models all nonradial modes as- 
sociated with the higher frequency maximum have mostly gravity 
character; they behave like p-modes only in the very outer enve- 
lope. 

This behaviour of 77 in the post-MS model can be explained by 
the effect of mode trapping. If there is a node of the radial displace- 
ment eigenfunction in the vicinity of the /j-gradient zone, a cavity is 
divided into two parts and a mode is trapped in the radiative helium 
core and in the envelope. In the envelope a mode behaves like a 



p-mode and energy trapped there is sufficient for the opacity mech- 
anism to be effective. In such a case, we have the positive value of 
77. The other situation is when there is no node in the vicinity of 
the /j-gradient zone. In this case a mode enters the radiative core 
with much higher amplitude and is very effectively damped. Then, 
a mode is stable. Therefore, two very close frequency modes can 
have opposite energetic properties; the one is unstable whereas the 
other is stable. 

We performed the instability analysis for models with masses 
from 7 M Q up to 20 M Q , the same as considered by Saio et al. 
(2006). Pulsational instability was found in post-MS models with 
masses greater than about 9 M Q up to 20 M Q . We did not consider 
models with masses higher than 20 M because it would require 
taking into account the effects of mass loss. As we did not consider 
models in the phase of the central helium burning, we evolved mod- 
els until the temperature in the center of a star reaches 10 8 K, i.e., 
before the onset of the 3a reactions. 

In Fig. 4, we show the instability domains on the HR diagrams 
marked as thick lines. Models with unstable radial and nonradial 
modes are depicted separately in the upper and lower panels, re- 
spectively. The left panels correspond to computations performed 
with the OP opacity tables and the right panels to those obtained 
with the OPAL data. The main difference is that with the OP data 
unstable modes appear from the very beginning on main sequence, 
i.e., from ZAMS for all masses but M = 20M o . With the OPAL 
tables we got instability from ZAMS for M < 1 1M and for higher 
masses it is shifted to the later phases on main sequence. As for 
pulsational instability of post-MS models, the difference between 
computations with the OP and OPAL data is subtle; with the OPAL 
tables the instability domain extends to slightly lower effective tem- 
peratures. 

Our results differ somewhat from those obtained by Saio et 
al. (2006) who adopted the OPAL data. Unlike these authors we 
did not have gaps in instability for the post-MS models with M = 
11 - 13M . Moreover, our instabilities do not extend to such low 
effective temperature as about 14 000 K. This may be due to dif- 
ferent chemical mixture adopted and/or differences in evolutionary 
and pulsational codes. 



3 PROPERTIES OF THE SPBSG MODES 

The asymptotic relations for high-order g-modes predict that at 
a given mode degree, C, periods of consecutive radial orders, n, 
should be equidistant (Tassoul 1980). It is clearly visible in the mid- 
dle panel of Fig. 5 where the value of AP, defined as 



AP = P„ 



P„-i 



(3) 



is plotted as a function of the frequency, v, for t = 1 and 2. We can 
see that the period spacings are generally constant and amounts to 
AP a 0.021 [d] for t = 1 and 0.012 [d] for 6 = 2, but there are local 
minima occurring when the instability parameter reaches maxima 
(upper panel of Fig 5). This behaviour of AP and its correlation 
with 77 is a consequence of mode trapping. 

It is also interesting to note that the period spacing between 
consecutive maxima of t/(v) is not constant and amounts to AP max = 
0.45 - 0.57 [d] for i = 1 and AP raax = 0.25 - 0.42 for I = 2 (the 
lower panel of Fig. 5). These maxima correspond to eigenmodes of 
the outer g-cavity. Therefore, in the case of oscillation spectra of 
the SPBsg pulsators one should look rather for irregularities corre- 
sponding to AP max . 

The phenomenon of mode trapping can be better explained by 
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Figure 6. The kinetic energy density (top panels) and differential work integral (bottom panels) of unstable (left panels) and stable (right panels) t ■ 
with close frequencies. The values of v and r\ are listed in panels. The same OP model as in Fig. 1 was used. 
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Figure 5. Correlation between the instability parameter, r] (top panel), pe- 
riod spacing, AP (middle panel) and AP max (bottom panel) in the frequency 
range of high-order g-modes with t = \ and 2 in the post-MS OP model 
with M=16 Me, log T eff = 4.340 and log L/L G = 4.730. 



the distribution of the kinetic energy density, E k , of modes and its 
radial and horizontal components. We adopted the usual normaliza- 
tion J E k dx = 1, where x = r/R. In Fig. 6, we plotted the kinetic 
energy density (upper panels) and differential work integral (lower 



panels) for two very close frequency modes with ( = 1, one of 
which is unstable (left panels) and the other stable (right panels). 
The unstable mode has a frequency of 0.172 c/d and corresponds 
to the local maximum of j], whereas the stable mode, with a fre- 
quency of 0.166 c/d, corresponds to the local minimum of tj (cf. 
Fig. 5). We considered the same OP model as presented in Fig. 1. 
As we can see, the kinetic energy density is higher for the stable 
mode than for unstable one. For both modes the kinetic energy den- 
sity is almost entirely confined to the radiative core where damping 
is extremely effective, what can be estimated from the values of 
dW/d\og T shown in the lower panels of Fig. 6. The difference is 
that the unstable mode has some addition of E k in the outer layers. 
It is more visible when one decomposes the kinetic energy density 
into the radial and horizontal components. 

Fig. 7 shows the radial (upper panels) and horizontal (bottom 
panels) components of E k for the same unstable (left panels) and 
stable (right panels) high-order g-modes as considered in Fig. 6. 
The horizontal component, which is related to g-mode behaviour, is 
much higher than the radial component for all high-order g-modes 
in SPBsg stars but the small radial component is crucial for insta- 
bility. The radial component is related to the p-mode behaviour and 
it is contained mainly in the envelope of the star. Modes perfectly 
trapped in both cavities have smaller total kinetic energy, hence the 
damping in the core is smaller and the higher radial component of 
energy in the envelope makes the opacity driving mechanism effec- 
tive. On the contrary, partially trapped modes have very high energy 
in the core and not sufficient energy in the envelope in order for the 
K-mechanism to be effective. 

It is important to add, that partial reflection of some oscillation 
modes always occurs in those specific stellar layers where the struc- 



© 2002 RAS, MNRAS 000,[Tj-?? 



6 J. Daszynska-Daszkiewicz, J- Ostrowski, A. A. Pamyatnykh 



i = 1 

freq = 0.166 c/d 
eta = -0.791 
logT eff = 4.340 
logL/L s = 4.730 




0.00 



7.5 7.0 6.5 6.0 5.5 5.0 4.5 7.5 7.0 6.5 6.0 5.5 5.0 4.5 

logT 



Figure 7. The radial (upper panels) and horizontal (bottom panels) components of the kinetic energy density for modes considered in Fig. 5. 



ture is characterized by sharp spatial changes. For gravity modes, 
the behavior of the Brunt-Vaisala' frequency is critically important. 
In our case, the reflection of some g-modes occurs at the top of the 
hydrogen gradient zone where the Brunt-Vaisala frequency has a 
deep minimum and where actual temperature gradient is changed 
from the adiabatic value to the radiative one (logT ss 7.14 and 
r/R » 0.144). Modes with small amplitudes (nodes) at this layer 
are reflected here. In other words, a mode trapping occurs in two 
regions, i.e., below and above of this layer. Due to the very high 
value of N in deep interior, the oscillation spectrum is very dense, 
but due to quite small values of N in outer part only few modes 
from such a dense spectrum (most of waves are confined in deep 
interior) have nodes close to the layer under consideration and are 
definitively trapped in two regions. The trapping in the outer part is 
critical for the instability of such modes. 

This is shown in Fig. 8 which demonstrates the reflection (or 
trapping) of some modes at the layer where the Brunt-Vaisala fre- 
quency has deep minimum. Again, two dipole gravity modes with 
close frequencies are considered. The trapped mode g282 with the 
frequency 0.171 1 c/d is unstable and the untrapped mode g 2 yo with 
the frequency 0.1794 c/d is stable. From the top panel we can see 
that 17 % of the kinetic energy of the trapped mode is confined 
in the outer part (above minimum of A'), whereas for the untrapped 
mode almost the whole kinetic energy (more than 99 %) is confined 



in the deep interior where this mode is strongly damped (see Fig. 6 
where stellar regions of damping and excitation are shown). Note 
that with standard normalization used in linear computations (the 
amplitude of the radial displacement on the surface is assumed to be 
1) the total kinetic energy of the trapped mode g 2 82 is 22 times less 
that those of untrapped mode g27o, therefore the gigi mode can be 
excited much easier. Lower panels illustrate the reflection (or mode 
trapping) phenomenon. The trapped mode is evanescent at the re- 
flecting layer (which is illustrated in the small panel) and has high 
oscillation amplitudes above this layer whereas the untrapped mode 
is almost totally confined in the deep interior and has not evanes- 
cent behavior at the minimum of the Brunt-Vaisala frequency. 



4 PROSPECT FOR MODE IDENTIFICATION 

As we could see, a number of modes can be excited in the B-type 
supergiant models. Due to the specific behaviour of the instability 
parameter, r], (cf. Fig. 5) one should rather not expect a constant 
period spacing in oscillation spectra of this new emerged type of 
pulsators. Therefore, let us now consider a possibility of identifica- 
tion of the SPBsg modes from multicolour photometry. 

If all effects of rotation on pulsation can be ignored, the com- 
plex amplitude of the light variations in the passband A for a star 
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Figure 8. The normalized kinetic energy (top panel) and the kinetic energy density (lower panels) of two dipole modes with close frequencies. The trapped 
mode g282 with the frequency 0.1711 c/d is unstable, and the untrapped mode g27o with the frequency 0.1794 c/d is stable. The partial reflection (trapping) of 
g282 mode occurs at the top of the hydrogen gradient zone (at log r/R = -0.840, r/R = 0. 144, log T = 7. 142), where the Brunt- Vaisala frequency has a deep 
minimum (dotted line, in arbitrary units). The position of this minimum is marked by the vertical line. For the trapped mode, 17 % of the kinetic energy is 
confined in the outer part above the reflection point. 



pulsating in a mode with the frequency to, the degree £ and az- 
imuthal order in is expressed as (e.g., Daszyhska-Daszkiewicz et 
al. 2002): 

SH A (i) = -1.0&6sY';'(i,0)b A ( (Di J + D %{ + D{ ( ) (1) 

where 



D 2 , = (2+r>(i-t), 



K, 



2 + 



w 2 « 3 \aiog(3^|^i) 



GM ) aiogg eff 



(2d) 
(2b) 
(2c) 



Here, e is the intrinsic mode amplitude, i is the inclination angle 
and G, M, R have their usual meanings. The term Df e describes the 
temperature effects, D 2 j - geometrical effects, and Dj . - the influ- 
ence of the pressure changes. The terms D'\ . and D^ , include the 
perturbation of the limb-darkening. The disc-averaging factor, b\, is 
the integral of the limb-darkening weighted by the Legendre poly- 
nomial. Derivatives of the monochromatic flux, ^(r c ff,logg), are 



calculated usually from static, plane-parallel atmosphere models. In 
general, these derivatives depend also on the metallicity parameter 
[m/H] and microturbulent velocity £,. Here, we rely on the NLTE 
atmosphere models of Lanz & Hubeny (2007). We use the fluxes in 
photometric passbands and nonlinear limb-darkening coefficients 
computed by Daszyhska-Daszkiewicz & Szewczuk (201 1). 

The / parameter describes the ratio of the bolometric flux per- 
turbation to the radial displacement at the level of the photosphere: 



^bol 



Re{sfY';'(9,ip)e- 



(3) 



The values of / is complex and it is obtained from linear computa- 
tions of stellar nonadiabatic pulsations. 

The positions of unstable modes with the mode degree, I, up 
to 3, on the diagnostic diagram employing the B, R Johnson filters 
are presented in Fig. 9. We chose these filters because they are close 
to those used on the BRITE-Constellation (BRIght-star Target Ex- 
plorer), the first forthcoming Austrian-Canadian-Polish space mis- 
sion which will perform two-colour photometry. Here, we show 
results obtained with the OP data. In the left panel of Fig. 9, we 
put unstable modes for the post-main sequence models with a mass 
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4>B " *R [ rad l 
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4>B " 4>R I rad ] 



Figure 9. The photometric diagnostic diagrams in the Johnson BR passbands with positions of unstable modes with I < 3 found in the post main sequence 
models. The left panel includes unstable modes of the post-MS models with amass of 16 M Q and the right panel - unstable modes of thepost-MS models with 
masses of 14, 15, 16, 18 and 20 M Q . The OP opacities, AGSS09 mixture and NLTE model atmospheres were adopted. 



of M = 16M Q and effective temperatures from log T e g = 4.384 to 
4.235. As we can see, separation of modes with different values of 
I is good. In the case of radial modes, the fundamental and first 
overtone modes follow different paths. This picture is not spoiled 
even if models in a wide range of masses are considered. In the 
right panel of Fig. 9, we put modes excited in models with masses 
M = 14, 15, 16, 18 and 20 M Q . The diagnostic photometric dia- 
grams for other pairs of passbands look quite similar but a combi- 
nation of the BR filters gives one of the best separation of modes. 

Using the OPAL data we get similar locations of modes with 
different mode degrees. The only difference is that there are less un- 
stable high-order g-modes modes because of the smaller maximum 
of r/(v) in the lower frequency for the OPAL pulsational models. 



5 CONCLUSIONS 

The goal of this paper was to analyse the pulsational instability of 
the B-type supergiant stellar models and to examine properties of 
their oscillation modes. We have showed that pulsational modes in 
the post-MS massive star models can be reflected at the chemical 
composition gradient zone developed above the radiative helium 
core and then the excitation in the Z-bump is effective. Both radial 
and nonradial modes are present. Our results are in contradiction 
with the results obtained by Saio et al. (2006) and Godart et al. 
(2009) who found that a presence of an intermediate convective 
zone (ICZ) in the hydrogen-burning shell is a necessary condition 
to get unstable modes in these hot evolved stars. 

As in the case of main sequence pulsators, the instability pa- 
rameter, rj, of post-MS models has maxima both at low and higher 
frequencies. In the low frequency range, the values of r/ reach many 
local maxima and minima which correlate with the period spac- 
ing. The two modes with very close frequencies can differ in their 
kinetic energy and stability properties. The kinetic energy and its 
distribution inside the star is different for unstable and stable high- 
order g-modes. For instability, the radial component of Et is cru- 
cial. These properties of the SPBsg modes are consequences of the 
nature of mode trapping in two cavities of the star, i.e., in the radia- 
tive helium core and in outer envelope. 



We obtained that equidistant pattern in the oscillation spectra 
of the SPBsg pulsators is rather not expected. This fact prompted us 
to explore a possibility of identification of the SPBsg modes from 
the multicolour photometry. We showed that the SPBsg modes with 
different degree I are located in well separated regions on the pho- 
tometric diagnostic diagrams. In particular, a very good separation 
was obtained for a pair of the BR Johnson passbands which are 
quite similar to those applied on the BRITE satellites. Therefore, a 
great hope is pinned on this forthcoming space project. 
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